verde.folha = function(da, mens="peso", ext="etanol")

{if(mens=="area"){

if(ext=="etanol") {

## renomeando a planilha ##

colnames(da)=c("amostra","tratamento" ,"volume", "area", "a470" , "a649" , "a665") 
############################################################################
########## Extrao com alcool 96% (lichtentaler et al, 1983) ##############



#####pigmentos em micromols por  ml########
ca.a=(13.95*(da$a665))-(6.88*(da$a649))
cb.b=  (24.96*(da$a649))-(7.32*(da$a665))
carot.c=((1000*(da$a470))-(2.05*ca.a)-(114.8*cb.b))/245



####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(da$area))
cb=((cb.b*(da$volume))/(da$area))
carot=((carot.c*(da$volume))/(da$area))


##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,ca,cb,carot)

#########Analise exploratria dos dados################
yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)

return(data.cl)}

#############################################################
########## Acetona 80%  0.1 - 0.5(Wellburn, 1994)  ###########
 if(ext=="acetona80.1"){

## renomeando a planilha ##

colnames(da)=c("amostra","tratamento", "volume","area", "a470" , "a663" , "a646") 

####pigmentos em micromols por  ml########
ca.a=(12.25*(da$a663))-(2.79*(da$a646))
cb.b=  (21.5*(da$a646))-(5.1*(da$a663))
carot.c=((1000*(da$a470))-(1.82*ca.a)-(85.02*cb.b))/198


####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(da$area))
cb=((cb.b*(da$volume))/(da$area))
carot=((carot.c*(da$volume))/(da$area))


##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,ca,cb,carot)

#########Analise exploratria dos dados################

yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)

return(data.cl)
}
########## Acetona 80% 1-4 nm(lichtentaler et al, 1983; Wellburn, 1994) ###########
 if(ext=="acetona80.2")
{
## renomeando a planilha ##

colnames(da)=c("amostra","tratamento", "volume","area", "a470" , "a663" , "a646") 

####pigmentos em micromols por  ml########
ca.a=(12.21*(da$a663))-(2.81*(da$a646))
cb.b=  (20.13*(da$a646))-(5.03*(da$a663))
carot.c=((1000*(da$a470))-(3.27*ca.a)-(104*cb.b))/229


####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(da$area))
cb=((cb.b*(da$volume))/(da$area))
carot=((carot.c*(da$volume))/(da$area))


##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,ca,cb,carot)

#########Analise exploratria dos dados################

yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)
return(data.cl)
}
#############################################################
########## acetona 100% (lichtentaler et al, 1983) ###########
 if(ext=="acetona100")
{
## renomeando a planilha ##

colnames(da)=c("amostra","tratamento", "volume","area", "a470" , "a645" , "a662") 

#####pigmentos em micromols por  ml########
ca.a=(11.75*(da$a662))-(2.35*(da$a645))
cb.b=  (18.61*(da$a645))-(3.96*(da$a662))
carot.c=((1000*(da$a470))-(2.27*ca.a)-(81.4*cb.b))/227


####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(da$area))
cb=((cb.b*(da$volume))/(da$area))
carot=((carot.c*(da$volume))/(da$area))
##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,ca,cb,carot)

#########Analise exploratria dos dados################

yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)
return(data.cl)
}
#################################################################
########## DMSO -range 0.1 a 0.5 nm (Wellburn, 1994) ###########
 if(ext=="dmso1")
{
## renomeando a planilha ##

colnames(da)=c("amostra","tratamento", "volume","area", "a480" , "a649.1" , "a665.1")

#####pigmentos em micromols por  ml########
ca.a=(12.47*(da$a665.1))-(3.62*(da$a649.1))
cb.b=  (25.06*(da$a649.1))-(6.5*(da$a665.1))
carot.c=((1000*(da$a480))-(1.29*ca.a)-(53.78*cb.b))/220

####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(da$area))
cb=((cb.b*(da$volume))/(da$area))
carot=((carot.c*(da$volume))/(da$area))
##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,ca,cb,carot)

#########Analise exploratria dos dados################

yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)
return(data.cl)
}
##########################################################
########## DMSO -range 1 - 4nm (Wellburn, 1994) ###########
 if(ext=="dmso2")
{
## renomeando a planilha ##

colnames(da)=c("amostra","tratamento", "volume","area", "a480" , "a649" , "a665")

#####pigmentos em micromols por  ml########
ca.a=(12.19*(da$a665))-(3.45*(da$a649))
cb.b=  (21.99*(da$a649))-(5.32*(da$a665))
carot.c=((1000*(da$a480))-(2.14*ca.a)-(70.16*cb.b))/220

####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(da$area))
cb=((cb.b*(da$volume))/(da$area))
carot=((carot.c*(da$volume))/(da$area))
##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,ca,cb,carot)

#########Analise exploratria dos dados################

yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)
return(data.cl)}}

{if(mens=="peso"){

#########      #########      #############    ########        #############
############################################################################
########## Extrao com alcool 96% (lichtentaler et al, 1983) ##############

if(ext=="etanol") {

## renomeando a planilha ##

colnames(da)=c("amostra","tratamento" ,"peso.fresco" ,"peso.saturado", 
"peso.seco", "peso.fresco.cl", "volume", "a470" , "a649" , "a665") 


###cra###
cra=((da$peso.fresco-da$peso.seco)/(da$peso.saturado-da$peso.seco)*100)

######peso seco cl= ((peso.fresco.cl * peso.seco)/peso.fresco))



peso.seco.cl=((da$peso.fresco.cl*da$peso.seco)/(da$peso.fresco))


#####pigmentos em micromols por  ml########
ca.a=(13.95*(da$a665))-(6.88*(da$a649))
cb.b=  (24.96*(da$a649))-(7.32*(da$a665))
carot.c=((1000*(da$a470))-(2.05*ca.a)-(114.8*cb.b))/245

####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(peso.seco.cl))
cb=((cb.b*(da$volume))/(peso.seco.cl))
carot=((carot.c*(da$volume))/(peso.seco.cl))

##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,cra,ca,cb,carot)

#########Analise exploratria dos dados################
yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)

return(data.cl)
}

#############################################################
########## Acetona 80%  0.1 - 0.5(Wellburn, 1994)  ###########
 if(ext=="acetona80.1") {

## renomeando a planilha ##

colnames(da)=c("amostra","tratamento" ,"peso.fresco" ,"peso.saturado", 
"peso.seco", "peso.fresco.cl", "volume", "a470" , "a649" , "a665") 

###cra###
cra=((da$peso.fresco-da$peso.seco)/(da$peso.saturado-da$peso.seco)*100)

######peso seco cl= ((peso.fresco.cl * peso.seco)/peso.fresco))

peso.seco.cl=((da$peso.fresco.cl*da$peso.seco)/(da$peso.fresco))


#####pigmentos em micromols por  ml########
ca.a=(12.25*(da$a663))-(2.79*(da$a646))
cb.b=  (21.5*(da$a646))-(5.1*(da$a663))
carot.c=((1000*(da$a470))-(1.82*ca.a)-(85.02*cb.b))/198

####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(peso.seco.cl))
cb=((cb.b*(da$volume))/(peso.seco.cl))
carot=((carot.c*(da$volume))/(peso.seco.cl))

##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,cra,ca,cb,carot)

#########Analise exploratria dos dados################
yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)

return(data.cl)}

########## Acetona 80% 1-4 nm(lichtentaler et al, 1983; Wellburn, 1994) ###########
 if(ext=="acetona80.2") {

## renomeando a planilha ##

colnames(da)=c("amostra","tratamento" ,"peso.fresco" ,"peso.saturado", 
"peso.seco", "peso.fresco.cl", "volume", "a470" , "a646" , "a663") 


##cra##
cra=((da$peso.fresco-da$peso.seco)/(da$peso.saturado-da$peso.seco)*100)

######peso seco cl= ((peso.fresco.cl * peso.seco)/peso.fresco))

peso.seco.cl=((da$peso.fresco.cl*da$peso.seco)/(da$peso.fresco))


#####pigmentos em micromols por  ml########
ca.a=(12.21*(da$a663))-(2.81*(da$a646))
cb.b=  (20.13*(da$a646))-(5.03*(da$a663))
carot.c=((1000*(da$a470))-(3.27*ca.a)-(104*cb.b))/229

####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(peso.seco.cl))
cb=((cb.b*(da$volume))/(peso.seco.cl))
carot=((carot.c*(da$volume))/(peso.seco.cl))

##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,cra,ca,cb,carot)

#########Analise exploratria dos dados################
yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)

return(data.cl)}
#############################################################
########## acetona 100% (lichtentaler et al, 1983) ###########
 if(ext=="acetona100"){ 

## renomeando a planilha ##

colnames(da)=c("amostra","tratamento" ,"peso.fresco" ,"peso.saturado", 
"peso.seco", "peso.fresco.cl", "volume", "a470" , "a645" , "a662") 


##cra##
cra=((da$peso.fresco-da$peso.seco)/(da$peso.saturado-da$peso.seco)*100)

######peso seco cl= ((peso.fresco.cl * peso.seco)/peso.fresco))

peso.seco.cl=((da$peso.fresco.cl*da$peso.seco)/(da$peso.fresco))


#####pigmentos em micromols por  ml########
ca.a=(11.75*(da$a662))-(2.35*(da$a645))
cb.b=  (18.61*(da$a645))-(3.96*(da$a662))
carot.c=((1000*(da$a470))-(2.27*ca.a)-(81.4*cb.b))/227

####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(peso.seco.cl))
cb=((cb.b*(da$volume))/(peso.seco.cl))
carot=((carot.c*(da$volume))/(peso.seco.cl))

##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,cra,ca,cb,carot)

#########Analise exploratria dos dados################
yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)

return(data.cl)}
#################################################################
########## DMSO -range 0.1 a 0.5 nm (Wellburn, 1994) ###########
 if(ext=="dmso1"){

## renomeando a planilha ##

colnames(da)=c("amostra","tratamento" ,"peso.fresco" ,"peso.saturado", 
"peso.seco", "peso.fresco.cl", "volume", "a480" , "a649.1" , "a665.1")


##cra##
cra=((da$peso.fresco-da$peso.seco)/(da$peso.saturado-da$peso.seco)*100)

######peso seco cl= ((peso.fresco.cl * peso.seco)/peso.fresco))

peso.seco.cl=((da$peso.fresco.cl*da$peso.seco)/(da$peso.fresco))


#####pigmentos em micromols por  ml########
ca.a=(12.47*(da$a665.1))-(3.62*(da$a649.1))
cb.b=  (25.06*(da$a649.1))-(6.5*(da$a665.1))
carot.c=((1000*(da$a480))-(1.29*ca.a)-(53.78*cb.b))/220

####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(peso.seco.cl))
cb=((cb.b*(da$volume))/(peso.seco.cl))
carot=((carot.c*(da$volume))/(peso.seco.cl))

##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,cra,ca,cb,carot)

#########Analise exploratria dos dados################
yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)

return(data.cl)}

##########################################################
########## DMSO -range 1 - 4nm (Wellburn, 1994) ###########
if(ext=="dmso2"){

## renomeando a planilha ##

colnames(da)=c("amostra","tratamento" ,"peso.fresco" ,"peso.saturado", 
"peso.seco", "peso.fresco.cl", "volume", "a480" , "a649" , "a665")


##cra##
cra=((da$peso.fresco-da$peso.seco)/(da$peso.saturado-da$peso.seco)*100)

######peso seco cl= ((peso.fresco.cl * peso.seco)/peso.fresco))

peso.seco.cl=((da$peso.fresco.cl*da$peso.seco)/(da$peso.fresco))


#####pigmentos em micromols por  ml########
ca.a=(12.19*(da$a665))-(3.45*(da$a649))
cb.b=  (21.99*(da$a649))-(5.32*(da$a665))
carot.c=((1000*(da$a480))-(2.14*ca.a)-(70.16*cb.b))/220

####### pegmentos em micromol por grama##########

ca=((ca.a*(da$volume))/(peso.seco.cl))
cb=((cb.b*(da$volume))/(peso.seco.cl))
carot=((carot.c*(da$volume))/(peso.seco.cl))

##########Matriz de resultados ############
amostra=da$amostra
tratamento=da$tratamento
data.cl=data.frame(amostra,tratamento,cra,ca,cb,carot)

#########Analise exploratria dos dados################
yy=c(data.cl$ca, data.cl$cb,data.cl$carot)
y=c(min(yy),max(yy)+0.5)

par(mfrow=c(1,3))
boxplot(ca~tratamento, data=data.cl, col="gray" , ylab= "clorofila a (micromol/g)", xlab="", ylim= y)
boxplot(cb~tratamento, data=data.cl, col="gray" , ylab= "clorofila b (micromol/g)", xlab="Tratamento", ylim= y)
boxplot(carot~tratamento, data=data.cl, col="gray" , ylab= "carotenoides (micromol/g)", xlab="",ylim= y)

return(data.cl)}}}}


